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l> ■ ABSTRACT 

C*~) ' A statistical analysis of the space-time distribution of absorption-line sys- 

o : 

, terns (ALSs) observed in QSO spectra within the cosmological redshift in- 



o 

systems (jRvabinkov et al 



terval z = 0.0-3.7 is carried ou t on the base of our catalog of absorption 



2003). We confirm our previous conclusion that 



the z-distribution of absorbing matter contains non-uniform component dis- 

, playing a pattern of statistically significant alternating maxima (peaks) and 

<3 ■ 

minima (dips). Using the wavelet transformation we determine the positions 
r ^ ■ of the maxima and minima and estimate their statistical significance. The po- 

X 

5-h ' sitions of the maxima and minima of the z-distributions obtained for different 

celestial hemispheres turn out to be weakly sensitive to orientations of the 
hemispheres. The data reveal a regularity (quasi-periodicity) of the sequence 
of the peaks and dips with respect to some rescaling functions of z. The same 
periodicity was found for the one-dimensional correlation function calculated 
for the sample of the ALSs under investigation. We assume the existence of a 
regular structure in the distribution of absorption matter, which is not only 
spatial but also temporal in nature with characteristic time varying within 
the interval 150-650 Myr for the cosmological model applied. 
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Figure 1. z-distributions of absorption-line sys tems observed in QSO sp ectra within the redshift interval z = 0.0 - 3.7: 
(i) 847 systems are ta ken from the catalogue by ljunkkarinen et al.l ll lQ9lh and (ii) 2003 systems - from the catalogue by 
iRvabinkov et al.l l|2003h (see text). 

1 INTRODUCTION 



Kaminker. Ryabinkov fc Varshalovich 



2000, here- 



We continue our previous studies (e.g., 
after Paper I) of the space-time distribution of absorption-line systems (ALSs) imprinted 
in spectra of quasars (QSOs). Actually, ALSs contain basic information on the distribu- 
tion of matter between the observer and QSOs as well as on physical processes occurred 
in different epochs of the cosmological evolution. In the first stage (e.g., Paper I) we used 
as indicators of matter only data on the resonance absorption doublets of C IV and Mg II 
observed in QSO spectra at cosmological redshifts z = 0.2 — 3.2 paying special attention to 
dimini shing of possible selection effects. Then, we explored considerably larger number of 



ALS s (IRvabinkov. Kaminker fc Varshalovichll200ll ; hereafter Paper II) basing on the catalog 



by ljunkkarinen. Hewitt fc Burbidgd (1199 ll ) within the extended range z = 0.0 — 3.7. It was 
shown that the ^-distribution of ALSs displays a pattern of alternating maxima (peaks) and 
minima (dips) which are statistically significant against a smooth dependence (trend). It is 
essential that the positions of the peaks and dips turned out to be independent (within sta- 
tistical uncertainties) of observation directions. This suggested that the derived distribution 
of absorption matter is not only spatial but also temporal in nature. 

We present the results of an extended statist ical analysis performed by different methods 
on the base of our catalog of absorption systems iRvabinkov et al. 2003 ). Fig. [1] demonstrates 
two z-distributions N&s(z) of ALSs containing absorption lines of heavy elements within the 
redshift interval z = 0.0 - 3.7. One of them is obtained using 847 systems from the catalog 
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by lJunkkarinen et al.l (119911 ) and the second one is based on 2003 systems registered in 



the spectra of 6 6 1 QS Os (emission redshifts within z e =0.29 - 4.9) from the catalog by 



Ryabinkov et al. 



(120031 ). All ALSs under consideration comprise lines of heavy elements and 
may include up to 20 - 30 absorption lines detected in different regions of QSO spectra 
(predominantly within the interval 3000 - 8000 A). We exclude ALSs consisting only of 
neutral hydrogen lines as well as damped Lya absorption systems (DLA). All redshifts Zj 
registered in spectra of each QSO and fallen into the velocity interval 5v =500 km s _1 are 
treated as a single absorption system with an averaged redshift z a = (Y^Li z j) / n a, where 
n a is a number of redshifts within the interval. Both distributions in Fig. [1] are obtained 
using so-called sliding- average approach which represents a set of consecutive displacements 
of the averaging bin A z = 0.071 along z-axis with a step S z = 0.01. 

A comparison of the two ^-distributions (Fig. [TJ reveals similar patterns of the peaks 
and dips relatively smooth curves. The positions of majority of the peaks and dips remain 
the same after the extension of statistics and some of them become more significant. There 
are only a few exceptions concerned with splitting of some initially single peaks into double 
ones and an appearance of new peaks (see Table 1). 

In Papers I and II the statistical analysis of ALSs was performed using the averaging 
bin A z = 0.071. This value was chosen with employing the x 2 -criterion as a measure of the 
most significant deviations of the ^-distributions obtained for different values A z from the 
hypothesis of the uniform distribution (with use of the trend subtraction) . On this way the 
results of statistical considerations could be sensitive to artificial non-uniformities induced 
by a separation of the whole sample of the redshifts into bins. To exclude sensitivity of our 
results to the effects of the averaging bins we employ in this paper predominantly out-of-bin 
statistical technique. In particular, in Section 2 we apply a continuous wavelet transform to 
the whole sample of ALS redshifts to study in details the peak-and-dips sequence. 

In Section 3, we demonstrate isotropy of the ^-distribution, i.e., approximate indepen- 
dence the peak/dip positions of observation directions. In Section 4, we focus on a presumable 
regularity (quasi-periodicity) of the ^-distribution. In Section 5, we examine properties of 
the one-dimensional correlation function. In Section 6, a polemic on the presence of period- 
icities in the distribution of QSO redshifts is outlined and possible selection effects in our 
analysis are briefly discussed. Conclusions and a short discussion of possible interpretations 
of the results are considered in Section 7. 
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2 WAVELET TRANSFORM AND NON-UNIFORMITY OF THE 
Z-DISTRIBUTION 

In this section, we use the wavelet transform for the analysis of the consecution of all 
redshifts under consideration within the interval z = 0.0 - 3.7 without subdividing the z- 
points into certain statistical bins. The wavelet transform is a way to disclose a difference 
of the redshift distribution from a smooth dependence and to reveal a set of condensations 
and depressions (peaks and dips) beyond the Poisson fluctuations. Actually, the wavelet 
analysis is effective when a signal has no evident periodicity and contains different sets 
of wave number (frequency) components in different space (time) regions. This method 
appears to be appropriate for the study of the non-uniformities visible in Fig. [1] and allows 
us to estimate statistical significance of the peaks and dips in the z-distribution of ALSs. 
Detailed information on the wavelet analysis and its application s in physics and astrophysics 



may be found in numerous books, reviews, and papers (e.g., 



Chui 



1992 



Astaf 'eva 



Dremin. Ivanov Nechitai'lo 



2001 



1996 . 



and references therein). 

We apply so called "Mexican hat" wavelet: 

u' 
~2 



*pK x ( x ,y) = A t 1 -u 2 ] ex P 



u 



u(x,y) = {x- y)K x 



(1) 



where x is the current one-dimensional coordinate, y is the translation of the wavelet along 
the coordinate axis, K x is an analog of the wave number for one-dimensional Fourier trans- 
form, and A = A(K X ) is a normalizing function. Usually, the value 1/K X is called a dilation 
or a scale parameter. The wavelet (OQ) is proportional to the second derivative of the Gaus- 
sian distribution exp[— u 2 /2]/\/2n taken with sign (— ). The zero-momentum of Eq. (TjQ) 
J dx ipK x ( x, v) = J dy ipK x {x, y) is equal to zero in accordance with the admissibility condi- 



tion (e.g., 



Chui 



1992|). 



Let we have a signal F(x), where F(x) G L 2 -space of normalized functions determined 
over the real axis x. Then, the coefficients of the continuous wavelet transform are 



S(K x ,y) = j dxF(x)^ Kx (x,y), (2) 

where the integration is carried out over the real axis. We obtain two-dimensional scanning 
of the signal, the wave number K x and the coordinate y being independent variables. 

In this paper, we obtain the coefficients ([2]) in an unconventional manner when the 
signal F(z) is treated as the consecution of the discr ete redshifts Zj (i = 1,2, ...) inside the 
considered interval (similar approach was suggested by lSlezak et al.lll993l ). We approximately 



present F(z) as the sum of normalized delta functions localizing positions Z{ of all points in 
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the sample: 



N 



F( z ) « £ ^ *(* - *); I J* \m\ 2 « N tot , (3) 

where ICf ^> 2ir is an upper boundary of all wave numbers [K z ^ KJf) involved in the 
analysis, N to t - is a number of all absorption systems in the sample, in our case N to t = 2003. 

Following Eqs. (T5]) and the wavelet transform of the function F(z) yields the sum 
over all redshift s: 

,2" 

(4) 



o Ntot 
Z7T 



S(K Z , z) = — z) = A(K Z ) J2 [1 " «?] exp 

2 i=i i=i 



2 



where = z) = (z, — z)i^ z , the variables K x and y in Eq. (|2J) are replaced by K z 
and z, respectively; in the second equality of Eq. (j4j) the factor of 2n /fCf is included in the 
normalizing function A(K Z ). 

Under the assumption of the uniform distribution of points Zi over the interval z = 0.0 - 
3.7 the value S(K Z , z) represents a random function obeying approximately to the Gaussian 
distribution (N to t 3> 1)- The mean value of the distribution /jl[S(K z , z)] is equal to zero: 

fi[S(K z ,z)] = - — £ / dziip Kz (zi, z) = 0, (5) 

where the upper limit is Z max = 3.7 and the mean squared deviation a(S(K z , z)) may be 
put equal to unity: 

a 2 [S(K z ,z)} = -±-£ / d ^Vl(^, *) = I- (6) 

"max JO 

The latter condition allows one to determine A = A(K Z ) (in general case A = A(K Z , z)). Us- 
ing Eqs. PJ, (jSJ) and neglecting terms of the orders of exp(— K 2 z 2 ) <C 1 and exp[— K^(Z max — 
z) 2 } <Cl one can approximately write 



A(K Z ) « r K * Z ™* . (7) 
1 2j y 3v^Nt ot 1 ; 

Such a normalization allows to measure naturally two-dimensional wavelet coefficients S(K Z , z 

in units of the standard deviation. 

To decontaminate the wavelet coefficients S(K Z , z) from noise (e.g. 



Romeo et al. 



2003, 



20041 ) at small scales (high K z ) and from the trend at large scales (low K z ) we exploit so- 
called spectral filtration. For this aim, one can calculate the energy spectrum of the wavelet 
transform 

E S (K Z )~ Y,S\K Zl z n ), (8) 
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where index n = 1,2, ...363 runs over the chosen interval (0.04 ^ z ^ 3.67) with a small step 
(S z = 0.01). 

If we admit that the distribution of S(K z ,z) is Gaussian (with zero mean value and 
unit dispersion) and all S(K z ,z n ) in Eq. (JSj) are statistically independent (see below), then 
the value Es(K z ) obeys to the ^-distribution with 363 degrees of freedom. In this case we 
can fix a certain tabulated value Xi-p &t a given significance level p (e.g., 3a, p = 0.0025) 
and select the region of K z ^ 28, where Es(K z ) exceeds this level (critical interval for the 
hypothesis of the random distribution of the redshifts). On the other hand, at K z < 14 
the value E S (K Z ) displays a sharp rise caused by the smooth dependence (trend) of the 
^-distribution (see Fig. [TJ. As a result we choose a representative interval 15 ^ K z ^ 28, 
which is consistent with significant variations of the z-distribution, and cancel S(K Z ) for all 
others K z . 

The upper panel of Fig. [2] illustrates the two-dimensional wavelet coefficients S(K z ,z) 
in the chosen interval of K z . One can see almost vertical alternating black and white strips 
corresponding to positive (enhanced) and negative (reduced) values of the wavelet coeffi- 
cients, respectively. Centers of the strips do not change for different K z and these centers 
are consistent with the centers of peaks z max and dips z m i n obtained earlier in Papers I and 
II (see Table 1). 

The lower panel of Fig. [2]illustrates how wavelet coefficients allow one to localize positions 
of the centers of the bands and estimate significance levels of the peaks and dips. The 
panel shows the z-dependence of the extremal wavelet coefficients S ext (K*, z) = S(K*,z), 
where the values K* = K*(z) are obtained numerically for each fixed z from the condition 
\S ext (K*, z)\ = max(\ S(K Z , z)\). The extremal wavelet coefficients represent the variations 
of S(K Z , z) in the most pronounced formj^. 

Positions of peaks z max and dips z min are determined as weighted mean values (centers 
of gravity) of the points S cxt (K*, z) forming the peaks or dips with respect to zero level. 
The accuracy of such determination is ±0.02. As it is seen from Table 1 majority of the 
values z max and z m i n calculated for the wavelet coefficients are close to the values z mSuX and 
z min obtained in our previous statistical studies. There are only a few exceptions (indicated 
in Table 1 by bold font) concerned with two cases of splitting of peaks into double ones at 



Note that, the results are not too sensitive to the form of presentation. Thus, fixing a certain value K z within the chosen 
interval (e.g., K z = 20) one can obtain the curve S(20, z) which is weakly different from the curve shown in Fig. [2] 
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Figure 2. Results of the wavelet transformation. Upper panel: two-dimensional coefficients of the wavelet transform S(K Z , z)\ 
black strips illustrate positive values (S > 0), white strips - negative values (S < 0). Lower panel: extremal values of S(K z ,z) 
satisfying to the condition \S e xt(K*, z)\ = max (\S(K Z , z)\), where K z , and K* belong to the intervals 15 ^ (K z , K z ) ^ 28. 



2 max = 0.79 and 3.16, and three cases of appearance of new peaks at z max = 2.31, 2.66 and 
3.41. The latter peaks are located on the decreasing part of the N\g(z) distribution (Fig. [1]) 
where statistics is essentially poorer. Note that the peak at z max = 2.66 ± 0.02 was obtained 
earlier by the statistical analysis of C IV doublets (first column in Table 1). 

It follows from Table 1 that, in general, there is a good agreement between the new and 
previous results. In addition, the significance levels of the peaks and dips estimated using 
the wavelet coefficients S ext (K*, z) are systematically higher than similar significance levels 
obtained in Papers I and II. 

A set of our numerical simulations of the wavelet coefficients S(K Z , z) produced for 
Poisson distribution of JVt t points Zi using Eqs. (j3j) and ((Tj) has shown that the mean 



8 A. I. Ryabinkov, A. D. Kaminker, D. A. Varshalovich 



Table 1. Rcdshift values z max (maxima) and z m in (minima) of the z-distribution of ALSs a \ 



299 C IV & 216 Mg II 

Paper I 
A r As(z) - bin statistics 




847 ALSs 
Junkkarinen et al. (1991) 
N\s(z) - bin statistics 




2003 ALSs 
Ryabinkov et al. (2003) 
Wavelet analysis 




Zmax 


Zmin 


Significance 
level 


Zmax Significance 
level 


Zmin 


Significance 
level 




;nificancc 
level 




0.10 


2(7 


0.22 2a 


0.11 


2(7 


0.22 


2(7 


0.44 


0.31 


3<T 


0.45 3a 


0.30 


2(7 


0.42 


3(7 


0.77 


0.58 


3(7 


0.79 3a 


U.DD 
0.81 


3a 
la 


0.68 
0.87 


1.0(7 
1(7 


1.04 


0.96 


2(7 


1.10 la 


0.97 


2(7 


1.05 


3(7 


1.30 


1.18 


1(7 


1.28 la 


1.16 


2.5a 


1.28 


1(7 


1.46 


1.35 


1(7 


1.44 2a 


1.38 


la 


1.46 


.3(7 


1.63 


1.54 


2a 


1.63 2(7 


1.55 


2.5(7 


1.62 


2(7 


1.78 


1.71 


la 


1.75 Ict 


1.68 


2(7 


1.77 


1(7 


1.98 


1.84 




1 Oft O^- 

i.yt> Zit 


1.87 


2.(7 


1 (17 

i.y / 


2.5(7 


2.14 


2.07 


la 


2.13 la 


2.07 


2(7 


2.15 


2(7 


2.45 


2.24 


2a 


2.42 3(7 


2.23 
2.38 


1.5(7 
1(7 


2.31 

2.49 


1(7 
1.5(7 


2.64 








2.61 


2(7 


2.66 


1.5(7 


2.86 


2.72 


3<T 


2.87 3a 


2.74 


2(7 


2.87 


2(7 




3.05 


2a 


3.16 2cr 


3.00 
3.16 


1.5(7 
1.5(7 


3.09 
3.22 


1.5(7 
1.5(7 




3.40 


3<T 


3.56 2ct 


3.32 
3.49 


1.5(7 
1.5(7 


3.41 

3.61 


1(7 

2.5(7 



a ) Bold values 2 max mark divergences between our previous and present results 



squared deviation a, or the variance a 2 , calculated for the wavelet coefficients turned out 
to be lower than unity. Actually, the random wavelet coefficients obtained at different z 
and K z are not independent. In particular, we have a 2 [S(K z , z)] = 0.56 at K z = 15 and 
a 2 [S(K z , z)] = 0.66 at K z = 28. Such a decrease of the variance results in an increase of 
the significance levels of S(K Z , z). Consequently, the significance levels of peaks and dips 
indicated in Table 1 may be regarded as minimal estimations. Additional calculations have 
shown that statistical dependence of S(K Z , z) at different K z and z in Eq. ([8]) leads to 
lowering of the boundary xl- p an d as a result to expansion of the critical interval of K z into 
the range of higher K z . Thus the chosen interval 15 ^ K z ^ 28 may be treated as minimal 
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Figure 3. Statistical comparison of wavelet transformations produced for the z-distribution of ALSs (thick lines), the Gaussian 
distribution of 2000 random z-points with periodical variations of the mean value (thin lines - 'Gauss'), the lognormal distri- 
bution simulated on the base of the same Gaussian model (thin lines — 'Lognormal'; see text). Upper panel: skewness, middle 
panel: excess (kurtosis), lower panel: probability of i-deviations in a dependence on K z within the interval 15 $C K z ^ 70. 
Dashed lines indicate the boundaries of critical regions at the significance level 3o\ 



critical interval or minimal representative interval appropriate to the alternative hypothesis 
of a regular origin of the peaks and dips in Fig. [2j 

In this section and hereafter we use the hypothesis of Gaussian distribution of the absorb- 
ing matter as basic one for estimations of the statistical significance of the peaks and dips 
and their regularity. It follows from the results discussed that the simple Gaussian model 
for the distribution of ALSs may be rejected on the significance level > 2a. However, it is 
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worthwhile to verify consistency of the ALS distribution with a wider class of Gaussian-like 
and non-Gaussian distributions. For this aim, we simulated two particular stochastic models 
(basic hypotheses) of the absorbing matter distribution: 

(i) Gaussian binned distribution of 2000 redshift points Zi with a periodical expectation /j,(z) 
and the variance equal to unity: 

P [x(z)\ = -L exp {- [x(z) - f,(z)) 2 12} , (9) 

V Z7T 



where x( z) = [n(z)—n ]/y/no, n(z) is a number of random points within a bin z±A z /2, A z = 
0.071, n is a mean number of points averaged over all bins; y,{z) = 2.5cos(27rz 16/Z max ), 
numbers 2.5 and 16 are chosen as an amplitude (significance level) and harmonic number 
of the variations; (ii) so-called lognormal random field suggested in literature as a plausi 
ble m odel for density fluctuations in the Universe (e.g. 



Coles fc Jones 



1991 



Bi&Davidsen 



19971 ). The lognormal distribution of the same number of points has been produced numeri- 
cally from the distribution (jHJ) by a transformation of the variable x(z) into y(z) = exp[x(z)]. 
Accordingly, the wavelet transforms S(K Z , z) have been calculated for both distributions us- 
ing simple modification of Eq. (J4]) adjusted to binned statistics. 

Fig. [3] displays a comparison of statistical symmetry properties of the wavelet transfor- 
mations produced for two models (i) and (ii) with S(K Z , z) calculated for ALSs. Three panels 
in the Fig. [3] represent (from top to bottom) the skewness, kurtosis (excess), and probability 
of sign deviations (criterion of signs) calculated in standard technique. One can notice that 
the symmetry properties of the wavelet coefficients obtained for ALSs and 'Gauss' model are 
statistically close. For both samples the skewness and excess are situated near zero for all 
values K z and the probability of ± deviations variates near the number 0.5 corresponding to 
symmetrical distributions. On the other hand, the 'Lognormal' model notably differs from 
both symmetrical distributions by its asymmetry. 

These results and results of our additional calculations carried out for the standard 
Gaussian field and its transformation into the lognormal random field confirm that the 
Gaussian (or Gaussian-like) distribution may be chosen as a background for the statistical 
analysis of ALS. Let us remind as well that in our consideration we treat all ALSs detected 
along each line of sight and fallen into the averaging velocity interval Sv (predominantly 
5v = 500 km/s) as a single redshift. That smoothes away all nonlinear stochastic fluctuations 
of the ^-distribution and approaches them to Gaussian-like random field. 



The redshift distribution 11 




redshift z 

Figure 4. z-dependences of the wavelet coefficients S C xt(K*, z) (see text) calculated for 24 hemispheres in the equatorial {upper 
panel) and Galactic {lower panel) coordinate systems. Upper panel: the hemispheres are determined by the right ascensions 
a c — 90° ^ a ^ q c + 90°, a c stands for hemisphere centers, and declinations —90° SC <5 sC +90°; the vertical axis represents a c 
in degrees plotted with the step 15° (l h ). Lower panel: the hemispheres are determined by the Galactic longitudes l c — 90° SC 
I $J l c +90°, l c stands for hemisphere centers, and latitudes —90° SC b Si +90°; the vertical axis represents l c in degrees plotted 
with the step 15°. 

3 ISOTROPY OF THE DISTRIBUTIONS 

Let us test isotropy (or anisotropy) of the ALS distribution comparing the z-dependences 
of the wavelet coefficients S ext (K*,z) calculated for different hemispheres. The procedure 
is similar to that used in Papers I and II for verification of a dependence of Nas(z) (see 
Section 1) on orientations of hemispheres. 

Fig. H] shows two sets of redshift dependences of the value S ext {K*, z) plotted for 24 
hemispheres. Both panels in Fig. H] may be treated as two-dimensional distributions of the 
wavelet coefficients S ext . The vertical axis on the upper panel (equatorial CS) represents the 
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right ascension of the hemisphere centers a c (in degrees), the horizontal axes corresponds to 
the redshifts z. The black strips display positive (enhanced) values of the wavelet coefficients, 
Scxt > 0' an d the white strips display negative (reduced) values of the wavelet coefficients, 
S'ext < 0. One can see the set of continuous vertical black and white strips with centers at 
the same values of z max and z min as indicated in Table 1, i.e., z max = 0.42, 1.05, 1.28, 1.46, 
1.97, 2.87, 3.09 and z min = 0.11, 0.30, 0.55, 0.97, 1.16, 1.38, 1.87, 2.23, 2.74, 3.00, 3.32. These 
values of z max and z min are independent of a c within statistical errors. The rest black and 
white strips are more fragmentary. However, the strips comprise segments which are notably 
longer than 180°. Such segments exceed possible effects of enhanced concentration of ALSs 
(clumps) at certain z in separate directions and can not be reproduced by overlapping of 
hemispheres. 

The vertical axis on the lower panel represents the Galactic longitude of the hemisphere 
centers l c (in degrees). It is clear from comparison of the upper and lower panels that the 
positions of the vertical strips as well as their width do not change (within statistical uncer- 
tainties) if we use the Galactic coordinates instead of the equatorial ones. Consequently, the 
wavelet coefficients do not depend substantially on the orientation of hemispheres. We calcu- 
lated the correlation coefficients for the ^-distributions in 12 pairs of opposite (independent) 
hemispheres and found moderate positive correlation at the significance level > 2 — 3a, 
probably due to variations of peak-and-dip profiles with relatively stable centers of gravity. 

On the other hand, our calculations of the wavelet coefficients within separate quadrants 
(a c ±45°) instead of hemispheres (a c ±90°) display more irregular black and white domains 
around certain z m3jX and z m \ n . As a result the strips-like picture (as in Fig. HI) turns out to 
be more fragmentary. 

Thus our wavelet analysis of the ALS redshifts gives evidence in favour of (i) the existence 
of the statistically significant ( > 2a) peaks and dips in the z-distribution of the ALSs at 
certain z max and z min and (ii) approximate independence z max and z min obtained for different 
hemispheres of hemisphere orientations. 

4 REGULARITY (QUASI-PERIODICITY) OF THE PEAKS AND DIPS 

The positions of the maxima z max and minima z m i n presented in Table 1 are distributed 
nonuniformly. For instance, it may be verified numerically that the power spectrum cal- 
culated for the whole sequence of redshift points (see below) displays no significant peaks 
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in a wide region of z-scales (^-periods) including all intervals between neighbour maxima 
or minima in Table 1. Therefore, to reveal some regularity (e.g., quasi-periodicity) in the 
sequence of the peaks and dips we need to use so-called rescaling functions r — f(z) which 
satisfy, e.g., to periodicity condition: 

Cat 1 - C = /OCi 1 ) - /CCJ = const, (10) 
where m = 1, 2, ... - is the continuous numeration of the maxima and minima. 

In Papers I and II we tested the quasi-periodicity of the ALS distribution using several 
simple rescaling (trial) functions of z including the function log(l +z) which attracts special 
attention in literature in search of a periodicity of the QSO ^-distribution (see Section 6). In 
contrast to these papers, our new calculations based on the extended statistical sample have 
not revealed a significant (at level > 2a) periodicity with respect to the same trial functions. 
Therefore, we perform the power-spectrum analysis using some special trial function (unless 
the different functions are discussed) which appears to be more appropriate for testing of 
the quasi-periodicity of the ALS distribution: 

7j = r(zi) = sinh" 1 ^^) - ai) a 2 ], (11) 

where i = l,2,...N to t is a numeration of all ALSs (N to t = 2003); a\ and a 2 are constants; 
7](z) is the conformal time defined by 



rife 



1 



y/n M (l + Zf + (1 - Q M ~ O0(l + Z) 2 + til 



dz. 



(12) 



where we choose f2 M = 0.3 and Q A = 0.7. 

The spectral power of the distributi on is calculated without the use of statistical bins 



A z (e.g., 



Karlsson 



1977 



P(k) 



Arp et alJll990L Paper I): 

. 2 

27rfcr,' ' 



E 

i=l 



COS 



+ 



E 

i=l 



sin 



(13) 



where Tj is given by Eq. ffTTl) . f = r(z Ntot ) — t{z\) is the whole interval of values Tj, and k is 
a harmonic number. The periodicity yields the peak in the power spectrum V = max (P(k)) 
with a confidence probability 

P=[l-ex V (-V)], (14) 

where the confidence level is defined with respect to the hypothesis of the uniform distribu- 
tion of the values Tj. 

Results of the power-spectrum analysis are sensitive to the averaging velocity interval 
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Figure 5. Power spectra P(fc), k is the harmonic number, calculated with using the trial function ||TTJ, l|12j l for the complete 
sample of the values jsj, i = 1,2, ...Ntot (see text). Upper panel represents P(fc) denned by Eq. H13I I at the velocity interval 
Sv = 500km s" 1 of AS associations (along each line of sight) into single redshifts 2; (see Section 1). The main peak corresponds 
to k = 16. Lower panel plots the dependence of the highest-peak value P(fc = 16) on the interval Sv. The horizontal dash lines 
specify the confidence levels fi = 0.95 (2cr; only on the upper panel), /3 = 0.998 (3a), f3 = 0.999936 (4cr), (3 = 0.9999994 (5cr), 
respectively. 

5v of the ALS associations into separate redshifts (see Sect. [T]). The upper panel of Fig. 
[5] plots the results of the calculations based on Eq. (113]) . where <5t>=500 km s _1 and the 
constants in Eq. (ITT]) are ai=1.26, a 2 =3.37. One can see that in the represented interval of 
k the only appreciable peak P(k) appears for k — 16 at the significance level > 4a. This 
peak corresponds to a periodicity in units of the function ( TIT]) with the period At = f/k = 
0. 199 ± 0.001. Note that the most significant peak is very sensitive to variations of the values 
a\ and a 2 at fixed Qm and (or variations of Qm and £7a at fixed ai and a 2 ). The deviations 
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Figure 6. Same as in upper panel of Fig. [5] but for 1778 ALSs sampled with elimination of all groups of QSOs with their 
ALSs located within lxl deg 2 square regions on the sky. 



of ai and a 2 from indicated values exceeding ±0.01 and ±0.28, respectively, reduce the peak 
P(/c) to the significance level < 3a. 

The lower panel of Fig. shows the dependence of the main peak amplitude P(k = 16) 
on the velocity interval Sv. It is seen that the maximum of P(k = 16) increases smoothly 
with decreasing of Sv and reaches the significance level > 5a at Sv =50-70 km s _1 , being 
basically robust to variations of Sv. Let us note that such an averaging of redshifts along 
each line of sight (up to Sv = 10,000 km/s) smoothes down possible effects of longitudinal 
clumping of absorbing matter. 

In addition Fig. [6] demonstrates effects of plausible clumping of ALSs in fields perpendic- 
ular to lines of sight. We calculate the same power-spectrum as in Fig. 0but with elimination 
of all groups of QSO spectra if their coordinates turn out to locate within any lxl deg 2 
square region on the celestial sphere. Similar significant peak P(k) at k = 16 as in Fig. [5] 
indicates weakness of transversal clumping effects in our analysis. 
The rescaling 



Alam et al. 



unction given by Eqs. ffTTl) and ffl2|) is not unique. For example, following 
(120041 ) we may choose a model independent ansatz for the Hubble parameter 



H{z) = H y/Sl m (1 + zf + A + Ml + z) + A 2 (1 



z 



,2. 



Omo = 8ttG p /(3#o), ( 15 ) 

where H and p are the Hubble constant and matter density at the present epoch (z — 0). 
We calculate P(k) replacing r« and f in Eq. (Tl3!) by the conformal time rji and fj = i](z-^ tot ) — 
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Tj(zi), where the rescaling function i](z) is given by Eq. (|T2|) with the denominator replaced 
by H(z)/Hq from ( |T5|) . Varying Qmo and the coefficients A , Ai, A 2 , and keeping the equality 
A + A\ + A 2 = 1 — ^mo, we also obtain the significant peak P{k) of the power spectra at 
^mo — 0.3 and A = — 0.61, A\ = 3.11, A 2 = —1.8 on the level > 4cr (slightly lower than 
that in Fig. E). 

The other example of a rescaling function is the mean comoving number density of 
absorbers associat ed with filaments (the elements of so-calle d large scale structure; see Sect. 
[7]) introduced by iDemiahski. Doroshkevich &: Turchaninovl (120001 ) in their Eqs. (2.2) and 
(2.15) with the function B(z) given by their Eq. (2.3a). Employing this rescaling function 
we obtain the peak of P(k) at the same k — 16 but even more significant (> 5<r) than the 
peak in Fig. [51 However, additional powerful peaks of P(k) (less significant than the main 
peak) also appear for this trial function. In this case the periodicity becomes more complex 
and requires a special consideration. 

For comparison, using Eqs. (|T3l) . (TTTj) and (|T2l we have calculated additionally the power- 
spectrum P(k) of the distribution of emission-line redshifts z e detected in the spectra of 661 
original QSOs within the interval z = 0.29-4.9. We have not revealed significant peaks (with 
significance > 3a) for a wide region of periods around At ~ 0.2. Thus we conclude that 
there is no straightforward coupling between possible periodicities of the ^-distributions 
of ALSs and original QSOs. Note, however, that the sample of QSOs is not statistically 
representative. 

To sum up the results of this section, we conclude that the appearance of significant peaks 
in the power spectra (calculated with different trial functions) may be regarded as additional 
evidence for reality of the peak-depression sequence in the ^-distribution of ALSs and its 
nonrandom origin. At the same time this quasi-periodical distribution does not correlate 
with the distribution of the original QSOs. 



5 ONE-DIMENSIONAL CORRELATION FUNCTION 



We determine a two-point correlation function £ (5t) for the sample of ALSs in an unconven- 
tional way. Specifically, the variable 5t substitutes a comoving distance r between pairs of 
objects (e .g., galaxies or clusters of galaxies) in the standard two-point correlation function 
£(rH e.g., 



Peebles 



1993 



Landy &: Szalay 



1993 



Kerscher. Szapudi fc Szalay 



200 







Jones et al 



2004 . and references therein). In our consideration 5t is a measure of the separation of two 
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arbitrary ALSs numerated by i and j (i,j = 1, 2, ...N tot , i ^ j) in units of the trial function 
Sr^j = |Ti — Tj\ introduced by Eq. ( TTTT) . Accordingly, the intervals of the conformal time 



given by (TT21) and the redshifts are Srjij 
we can write: 



\rji — r]j\ and fejj 



A/"obs(^ 



respectively Thus 



(16) 



A4im(5r) 

A/"obs(^ r ) is the number of observed pairs of ALSs separated by the variables St^ within 
the range St ± A T /2, where A r = 0.07 is the chosen bin width. A4im( ( ^' r ) is the calculated 
number of cross pairs between the real sample of ASs and the points of a random sample 
simulated for the sa me intervals of z or t(z) and the same smoothed function (trend) as 
the real sample (e.g. jMo. Jing fc Bornerlll992l ). The both samples (real and simulated ones) 
have the same number of redshift points. 

Note that the variable St in Eq. ( fl6l) may be considered as a measure of the temporal 
separation of any two epochs in units of the scaling function (IH|) . Actually, in our case 
the temporal interpretation of £(5r) is more appropriate than the spatial one. Unlike the 
standard approach to the spatial two-point correlation functions (see below) the function 
(Fl6|) is based on the conception of a single reference center r = (z = 0) concerned with the 
observer. Accordingly, all sampled points (redshifts) Ti(zi) are distributed over concentric 
spherical layers (bins) independently of the directions to their original QSOs. We count up 
all redshifts Tj inside a concentric layer (r ± A T /2) with respect to any chosen point Tj (out 
of this layer). All redshifts within a layer are treated as equivalent despite of various spatial 
(comoving) distances between them. 

On the other hand, the function (fl6l) differs from the correlation functions calculated by 



counting only line-of-sight pairs of 



the o 



Quashnock. Vanden Berk fc York 



j jects with different mutua l comoving distances (e.g., 



1996 . 



Broadhurst fc Jaffe 



20001 ). although, the results of 



both approaches may turn out to be quite consistent. 

The upper panel of Fig. [7] displays the one-dimensional (two-point) co rrelation function 
ffT6l) in units of the appropriate Poisson error cr(£) calculated as follows (e.g., 
19911 ): 

£(5t) KUSt) - KUSt) 



Peacock fc Nicholson 



T 



where 



\M4im(<5 



(17) 



(18) 
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Figure 7. Upper panel: one-dimensional correlation function §(<5t)/(t(§) (see text) versus the interval 8r (0.1 < St < 1.9) 
between the components of pairs of ALS redshifts (sampled at Sv = 500km s _1 ); t(z) is the trial function given by Eq. 11H . 
Lower panel: power spectrum P^(fc) calculated for fi = fi(r) = £(<5t)/o"(£) according to Eq. (1191 : k is the harmonic number. 



One can see the sequence of positive and negative peaks with the significance > 4a 
with respect to zero level. Let us notice presence of the long range ordering in dependence 
of £(5r)/cr(£) on dr. We regard it as a consequence of the peak-and-dip structure in the 
z-distribution of ALSs (see Sect. [2]). 

The lower panel in Fig. [7| represents the result of the power spectrum analysis performed 
for the value \i = li(6t) = £(5r)/cr(£) according to the equation: 



Pm(*) 



1 



Nn 



m=l 



COS 



27rk 8t„ 
5f 



N u 



m=l 



sin 



5f 



(19) 
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where f m = /x m - < /x >, /x m = /x(5r m ), < /x > = N' 1 is the mean value of /x, 

D = < /i 2 > — < /x > 2 is the variance of /x, and < /x 2 > = A*"" 1 Xlm M Mm i s the mean 
squared value of /x; the values £r m run over a set of points m = 1,2, ...A^ of the variable 
#r from 0.1 to 1.86 with the step 0.01. Thus = 176 is the full number of the points and 
8t = 1.76 is the whole interval of Sr m variations under consideration. Note the strong single 
peak of P^(fc) at k ~ 9 on the significance level well exceeding 5a. This peak manifests the 
periodicity of the correlation function with the period At = 0.199 ± 0.001, which is in a 
good agreement with the periodicity found in Sect. HI 

The one-dimensional correlation functions £(5t) and appropriate power spectra ( fl9l) were 
calculated also for ALS samples with different values of Sv and Av, where 8v is the averaging 
velocity interval of ALS associations into single redshifts (see Sects. [T]andH]) and Av is the 
minimal velocity shift (along each line of sight) adopted for the sample of ALSs relative to 
the original QSO emission redshifts (z e ). In that way we exclude all ALSs with z a within 
the region associated with their own QSOs, i.e., at z e ^ z a ^ z e — Av (1 + z e )/c. The 
values 8v and Av were changed within the intervals 10 - 10,000 km/s and - 10,000 km/s, 
respectively. The results obtained are similar to that presented in Fig. [7J Thus the results 
are quite robust to the variations of ALS samples. 



For comparison, we have performed additional calculations o 



Landy fc Szalay 



1993 



Jones et al. 



: the standard two-point 



20041 ) . where r is the co- 



correlation function £(r) (e.g., 
moving distanc e between t he co mponents of ALS pairs. The values r were calculated using 
the formulae by iRoukemal (120011 ) with equatorial coordinates aj, and 5i appropriate to ALSs 
(Section |3J). For convenience, we employed the Hubble parameter H(z) given by Eq. ()15p 
with Q M0 = 0.3 and the same coefficients A , A\ and A 2 as indicated in Section HI Let 
us remind that the rescaling function ( [121) with the model dependence H(z) displays the 
significant periodicity of the ALS redshifts (see Section 4). We tested a periodicity of the 
function £(r) (in a similar way as for £(<5t)) using an analogue of Eq. ( fl9l) with the variable 
5r] = H r/c instead of 5r. As a result, we have not revealed any significant peaks of P At (fc) 
and, consequently, significant periodicities of the spatial correlation function £ (r) . 

Fig. El shows two sets of the values £(5t)/(t(£) as a function of 5t calculated for 24 
hemispheres, i.e., the step of consecutive rotations is also 15°. Similar to Fig. Hlthe upper 
and lower panels of Fig. [8] may be treated as two-dimensional distributions in the equatorial 
and Galactic coordinate systems, respectively. The black strips display the positive values 
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Figure 8. The correlation function £(<5r) of absorption systems calculated in 24 hemispheres for the equatorial (upper panel) 
and Galactic (lower panel) coordinate systems, respectively; £ > - black strips and £ < — white strips. Upper panel: the 
vertical axis indicates the right ascension of the hemisphere centers ct c (in degrees). Lower panel: the vertical axis indicates 
the Galactic longitudes of the hemisphere centers l c (in degrees). The horizontal axes show the interval of St between the 
components of ALS pairs. 

of the correlation function (£(5r)/cr(£) > 0) while the white strips display the negative ones. 
One can see again (cf. Fig. H|) that there is the set of continuous vertical black and white 
strips with centers at the same values 5t as on the upper panel of Fig. [7J It follows from 
comparison of the upper and lower panels that the positions of the vertical strips as well as 
their widths do not change essentially with transition from the equatorial to the Galactic 
coordinates, i.e., with changing of the axis of hemisphere rotations. Consequently, the two- 
point correlation function (similar to the wavelet coefficients in Sect. [3]) is weakly sensitive 
to orientations of the hemisphere. 
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6 DISTRIBUTION OF QUASARS AND SELECTION EFFECTS 



The first peak at z — 1.9 5 in the ^-distribution of emi ssion and absorption line redshifts in 
QSO spectra was found by lBurbidge fc Burbidgd (119671 ). Then the quasi-periodical sequences 
of pe aks and dips in the ^-distribut i on of quasar em i ssion l ines were marked by many authors 



Burbidgd 11968. 



Cowan 



1969 



Karlsson 



1971 



1977 



1990 



Fang et al. 



1982 



Arp et al. 



199CH ) . Furthermor e, several maxima and minima were detected in the distribution of absorp- 



tion systems (e.g., 



Chu & Zhu 


1989. 


Arp et al. 


1990) 



availability of the periodicities in the redshift distribution of QSOs or QSO-galaxy pairs have 



of the periodicitv does reallv exist (e.g., Karlsson 


1990. 


Khodvachikh 


1990. 


Arp et al 


1990. 


Arp. Fulton & Roscoe \ 


2005. 


Burbidge & Napier 


2001. 


Napier & Burbidge 


2003. Bu 


rbidge 


2003. 


Bell 


2004. 


and r 


eferences therein) and their opponents 




Box & Roeder 


1984. 


Basu 


1985 


2005 


. Scott 


1991 


. Hawkins. Maddox & Merrifieldl 


2002, 


Tang & Zhang 


2005). 



who argue that the periodicity and the very set of peaks and dips in the QSO distribution 
are results of various selection effects. 

At least a part of the authors standing for reality of the periodicity under discussion 
regards that as evidence for the hypothesis that QSOs (or essential part of them) repre 



sent s ome objects ejected from the nuclei of nearby active galaxies (e.g., I Arp et al. 



2005 



Burbidge &: Napier 



2001 



Bell 



199(1 . 



20041 . and references therein). The redshifts of QSOs 



being larger than the redshifts of their parent galaxies have been assumed to have a non- 
cosmological "intrinsic" origin and display a set of preferred (discrete) quasi-periodical val- 
ues. To our knowledge, there are two models discussed in literature which suggest different 
periodical sets of preferred redshifts. The first model has been applied mainly to QSOs and 



marginally to ALSs in QSO spectra (see above). That is desc ribed by so ca 



formula for the redshi 



Arp et al. 



1990 



2005 



t periodicity Alog(l + z f ) = 0.089 (e.g., 



Burbidge fc Napier 



was proposed by Bell (e.g., 



Bell 



2002 



2001 



Ka ' 



Napier fc Burbidge 



ssonl 



1971 



led Karlsson 



1977 



1990, 



20041 ) for the "intrinsic" redsh i 



20031 ). The second one 
ts of QSOs and ex- 



Bell fc Comeaul (120031 ). According to 



tended on a set of preferred redshifts of galaxies by 
both scenarios, the formation of ALSs on the lines of sight to the QSOs should have also 
occurred in the local vicinity of the Universe. 

To test a possible correlation of the ALS distribution with expectations of the first model 
we have calculated the power spectrum of the sample of ALSs with respect to the trial 



22 A. I. Ryabinkov, A. D. Kaminker, D. A. Varshalovich 

function log(l + z). We have found no significant power peaks which would have a chance to 
be related with the period 0.089 or with a multiple value of it. According to the second QSO 
ejection model, appearance of the preferable z max and z min in the ALS distribution could be, 
in principle, interpreted as effects of an additional set of discrete "intrinsic" redshifts referred 
to galaxies. To test this statement one can compare the values of z max from the last but one 
column in Table 1 with the set of intrinsic redshift components Zj o\N, m\ (where N an d 
m are some quantum numbers) defined for galaxies by Eq. (Bl) of iBell fc Comeaul (120031 ). 
Such a comparison shows that the set of the peaks z max given in Table 1 are not consistent 
with the intrinsic redshifts of galaxies. Only three values z max turn out to be equal (within 
the uncertainty ±0.02) to the values z i( j from the set of 30 numbers at N ^ 8 within the 
interval 0.1 ^ z ^ 3.7. Thus we have found no traces of consistency between our results and 
the hypotheses of non-cosmological "intrinsic" redshifts of QSOs. Accordingly, in contrast 
with the papers quoted above , we assume that the redshift of the absorption lines z a is 
cosmological and the observed ALSs are associated with ionized gas in intervening galaxies 
or clusters of galaxies at cosmological distances along lines of sight to original QSOs. 

The results of this paper may be biased by selection effects appropriate to the sample of 
ALSs. Taking it into account we undertook some special efforts in Paper I to minimize such 
effects in the statistical analysis. As it is shown in Table 1 the results obtained in Paper I 
are consistent with the results of the present paper. In Section [5] we discuss the effects of 
excluded regions of z (within Av) near QSO emission redshifts. Such exclusions allow us to 
avoid at least a part of selection effects related to so called associated regions of QSOs (e.g., 
Quashnock et al.lll996l ). 

Furthermore, we failed to find similar periodicity (at significance level > 3a) in the z- 
distribution of the 661 original QSOs as it is revealed in the distribution of ALSs registered 



in their spe ctra (see Section HI). 



QSOs (e.g., 



Box Roeder 



1984 



t separates possible selection effe cts biasing the statistics of 



Basu 



2005 . 



Tang &: Zhang)l2005l ) and those peculiar to the 



sample of ALSs involved in the analysis. It is worthwhile to note that the characteristic scale 
of redshifts z which we consider in this paper (~ 0.18) is essentially smaller than th e scale 



0.48) discussed for QSOs and QSO-galaxy p airs ( e .g.. iBurbidge fc Napier 



2001 or the 



^-period (~ 0.62) found for a sample of QSOs by iBelll (120041) . In part i cular , the narrowness 
of the peaks and dips obtained here makes the analysis of Basul ( 19791 . 1983 ) less relevant to 
our results. 

In Fig. [9] we test at least a part of possible selection effects concerned with a priori differ- 
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Figure 9. A-distribution (spectroscopic completeness) of resonance doublets and single absorption lines of 15 ions which 
contribute mostly in the distribution of ALSs at different z; bin size of the histogram is AA = lOOA. 



ent probability to register absorption lines in QSOs spectra at different values of wavelength 
A within the observational interval (3000 - 8000 A). We choose 15 single lines or resonance 
doublets of the ions (C II, C IV, Mg I, Mg II, Si II, Si III, Si IV, N V, Al II, Al III, Fe II) 
most representative in the sample of ALSs under consideration and determine a number 
of wavelengths A = Ao(l + z a ) observed for each of them, where Ao is a laboratory value. 
Fig. represents a histogram of the absorption lines as a result of count of the lines within 
bins AA = lOOA. One can see quite smooth histogram (so-called spectroscopic complete- 
ness) with two appreciable peaks between 4000 - 4200 A and 6300 - 6600 A (at significance 
level > 2a). These peaks are obviously not enough to explain all set of peaks in Fig. [1] 
as pure selection effects, i.e., as a consequence of peaks in the spectroscopic completeness. 
Although this conclusion does not completely exclude other selection effects in our statistical 
consideration. 



7 CONCLUSIONS AND DISCUSSION 

The main results of the statistical analysis of 2003 absorption-line systems (ALSs) in the 
redshift range z = 0.0 — 3.7 can be summarized as follows: 

(1) The z-distribution of ALSs displays the statistically significant pattern of alternating 
maxima (peaks) and minima (dips) relative to a smooth curve. The positions of the maxima 
and minima (centroids of the peaks and dips) in the ^-distribution are given in Table 1 
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Figure 10. Comoving distance intervals AR (left vertical axis) and appropriate time intervals AT (right vertical axis) 
between successive neighbour peaks in the 2-distribution of ALSs calculated as the intervals Arj multiplied by c/Hq and 1 /Ho, 
respectively (see text), and referred to appropriate redshifts z. The dotted lines illustrate ±<r deviations. 



with an accuracy ±0.02. These positions are predominantly not sensitive (within statistical 
uncertainties) to orientations of the hemispheres under statistical consideration. However, 
the correlation coefficients between the ^-distributions calculated for 12 pairs of opposite 
(independent) hemispheres turn out to be moderately significant ( > 2— 3a) as a consequence 
of variations of the peak-and-dip profiles at relatively stable weighted centers of the peaks 
and dips. 

(2) The sequence of peaks and dips reveals a certain regularity. The power spectrum 
calculated according to Eq. ({13]) with the rescaling function t(z) given by Eqs. (TTTj) and 
ffT2"j) displays the peak for the harmonic number k — 16 (period At = 0.199 ± 0.001) at 
the significance level exceeding 4cr (relatively to the hypothesis of the uniform distribution 
of ALSs over t(z)). The rescaling function t(z) may be treated as a measure of temporal 
separations of the epochs with different z. Still more prominent peak at the same period At 
arises in the power spectrum calculated for the two-point correlation function £(5t) with 
the use of Eqs. (tlB) and (tl9|. 

(3) The obtained distribution of ALSs is likely to be coupled with the appearance of 
alternating pronounced (peaks) and depressed (dips) epochs in the course of the cosmological 
evolution, i.e., with the existence of some (r elatively weak ) spatial-temporal wave process. 



According to the cosmological principle (e.g. iPeeblea Il993l ) similar wave-like process would 



be observed from any spatial-temporal points in the Universe. 
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(4) Fig. [10] demonstrates the intervals of the comoving distance, AR = cAtj/Hq, and the 
appropriate time, AT = Arj/Ho, between neighbour peaks in dependence on the current 
redshifts z, where rj(z) is given by Eq. (fl2l . At] runs over all successive intervals t](z l ^ x ) — 



v{4 



I = 1,2, ... - numerates the peaks in the ^-distribution of the ALSs (from top to 



bottom in Table 1)q It follows from Fig. [10] that the characteristic scales of AR and AT 
vary in the ranges AR = (50 - 200) h^ 1 Mpc and AT = (150 - 650) h^ 1 Myr, respectively, 
where h = H /100 km(s Mpc) -1 . 

Let us note that our treatment of the results obtained does not contradict to the exis- 



tence of the Large Sea 



Shandarin fc Zeldovich 



2004 



Jones et al 



e Structure (LSS) of the matter distribution in the Universe (e.g., 



1989 



Einasto et al. 2002 



2004, lEisenstein et al 



Tago et al 



2002 



Demiahski fc Doroshkevich 



20051 and references therein) . The spatial-temporal 



variations discussed here may be superimposed upon the process of formation and evolution 
of the LSS elements. It was noted in Sect. [3] that the strip-like picture represented in Fig. H] 
becomes more fragmentary if the ^-distributions are calculated for the quadrants (a c ± 45°) 
instead of the hemispheres. One can assume that more fragmentary patterns obtained for 
smaller sectors of observations appear not only as a result of a reduction of ALS statistics 
but also due to an interplay between LSS and the variations discussed here. In any case the 
possibility of such an interplay deserves a special investigation. 

Extended discussions of the effects of the periodicity in the redshift distribution of galax- 
ies on a s cale about 130 h~\ Mpc was initiated in literature by the results of a pencil-beam 



survey of 



Broadhurst et al 



(1l990l ). These results also admit both interpretations either as a 



pure spatial quasi-periodic pattern of LSS consti t uents (set o: 



Kaiser fc Peacock 



1991 



van de Weygaert 



1991 



Dekel et al. 



clumps or walls and yoids; e.g. , 



1992 



Yoshida et al. 



2001 



and 



references therein) or as a temporal sequence of pronounced and depressed epochs becoming 
apparent in the matter distribution. The latter treatment is closer to our interpretation of 
the results represented here. However, our statistical approach does principally differ from 
the pencil-beam consideration. 

The appearance of alternate pronounced and depressed epochs may b e associated wit h 



Liu&Hu 



19981 ). 



intrinsic quasi-periodical properties of galaxies and/or their clusters (e.g., 
On the other hand, it can be explained as an apparent pattern of density variations in 
terms of the cosmological scenarios considering scalar (scalar-tensor) fields and/or phase 



2 Note that the values AR and AT were calculated with normalization of dz in Eq. H12H to unity 
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transitions at different evolutionary stages. Such effects may lead to temporal oscillations 



of the Hubble parameter H(z) (e .g. 



Morikawa 



1990, 



1991 



Hill. Steinhardt fe Turner 



1994 



1994 



Busarello et al. 



1994 



Perivqlaropoulos fe Sourdis 120021) and/or effective g ravitational constant G 



2005) . peculiar velocity field ( 



Gonzalez et al. 



Hill et a 



Dodelson. Kaplinghat. Stewart 



2001 



Banerjee. Pavon fe Sen 



19911) . cosmological constant A (e.g.. 



1990, 



Sisterna fe Vucetich 



2003 



Davidson 



Damour fe Polyakov 



2000j), and, as a consequence, to apparent oscillations 



of the matter distribution with the redshift. 

The rapid development of ground-based and space-born observational facilities in the 
optical, infrared, and near-ultraviolet ranges in recent years gives confidence that the hy- 
pothesis of the existence of the alternate epochs with enhanced and reduced concentrations 
of absorbing (or luminous) matter in the Universe may be tested in the nearest future. 
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